Compressibility of an Ultracold Fermi Gas with Repulsive Interactions 
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Fermi gases with repulsive interactions are characterized by measuring their compressibility as a 
function of interaction strength. The compressibility is obtained from in-trap density distributions 
monitored by phase contrast imaging. For interaction parameters kpa > 0.25 fast decay of the gas 
prevents the observation of equilibrium profiles. For smaller interaction parameters, the results are 
adequately described by first-order perturbation theory. A novel phase contrast imaging method 
compensates for dispersive distortions of the images. 
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I. INTRODUCTION 



Experiments with ultracold atoms explore many-body 
physics with strong interactions. They have demon- 
strated long-predicted phenomena like the BEC-BCS 
crossover [l| and Lee-Huang- Yang corrections to the en- 
ergy of degenerate gases [2tl] • Experiments have also ex- 
plored novel quantumphases like fermions with unitarity 
limited interactions [J H|, population imbalanced Fermi 
gases [1, 0| and Hubbard models in optical lattices [1, ■ 
More recently, they have been used to provide precision 
tests of many-body theories [f|. Usually, interactions 
in ultracold gases are fully described by the scattering 
length, which is a zero-range approximation greatly sim- 
plifying the theoretical description. This approximation 
is valid since the diluteness of the atomic gases implies a 
particle separation much larger than the short range of 
the van-der-Waals interactions. This almost exact char- 
acterization of the interactions by a single parameter, the 
tunability of interaction strength, and precise experimen- 
tal control over cold atoms systems have made them an 
ideal testbed for many-body quantum calculations. 

A new level of quantitative comparison between the- 
ory and experiment was recently reached by careful mea- 
surements of density profiles from which the equation 
of state could be determined. These techniques were 
first proposed by Chevy [l(| and Bulgac [ll| and imple- 
mented by Shin 12 1 . Further improvements [1, [H, Il3l - |l6j | 
resulted in impressive accuracy without adjustable pa- 
rameters. These results hinge on accurate measurements 
of the equilibrium atomic density distribution. Since all 
cold atom systems are in a metastable phase, this requires 
a favorable ratio of lifetime to equilibration time. 



Long lifetimes and strong interactions were realized in 
Fermi gases with strong attractive interactions since the 
decay to lower-lying molecular states is suppressed by the 
Pauli exclusion principle [l7j . This is different for repul- 
sive interactions which are realized on the so-called upper 
branch of a Feshbach resonance where decay is always 
possible into the so-called lower branch which consists 



of weakly bound molecular states with binding energy 
H 2 /ma 2 with a being the scattering length. 

For bosons, the first beyond mean field correction, 
the so-called Lee-Huang- Yang term, could be observed, 
but corrections were necessary to account for the non- 
equilibrium profile since the time to sweep to strong in- 
teractions was not long compared to equilibration times 
and inverse loss rates Q. Here we study fermions with 
repulsive interactions. They have been the focus of much 
recent work due to the prediction of a phase transition to 
a ferromagnetic state for sufficiently strong interactions 
[l8l [l9j . Recent experimental [2(| HH and theoretical 
studies (22[ addressed the competition with strong decay 
to the lower molecular branch. 

As expected we find only a limited window for metasta- 
bility where we can observe equilibrated clouds and char- 
acterize the repulsive interactions by obtaining the com- 
pressibility from observed profiles. We observe the linear 
term in kpa corresponding to mean field energy for the 
first time in density profiles, hp is the Fermi wave vector. 
In contrast to a Bose-Einstein condensate, the mean field 
energy is smaller than the kinetic energy and also com- 
petes with thermal energy, and is therefore much more 
difficult to observe. The signal-to-noise ratio (and some 
heating) prevented us from discerning the second-order 
interaction term which is the Lee- Huang- Yang correction 
for fermions. 

Our work features one technical novelty, a novel im- 
plementation of phase-contrast imaging to address dis- 
persive distortions of the cloud. All studies mentioned 
above, with one exception [12], were conducted using res- 
onant absorption where dispersion (an index of refraction 
different from one) is absent. However, this severely lim- 
its the cloud size and number of atoms to optical densities 
of a few. Phase-contrast imaging has many advantages. 
It can be applied to clouds with much larger optical den- 
sities by adjusting the detuning. Due to coherent for- 
ward scattering, the heating effect per detected signal 
photon is reduced by potentially a large number (which 
is equal to the resonant optical density divided by four 
[HI]). This can e.g. be used for repeated nondestruc- 
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tive imaging. However, for precision studies of density 
profiles, small dispersive distortions of the density profile 
cannot be neglected. Previous work including Ref. [l2[ 
was not sensitive to this effect. We have developed an 
experimental technique to correct for dispersion. 



II. EXPERIMENTAL SETUP 

A spin polarized Fermi gas of 6 Li in the 
\F = 3/2; rap = 3/2) state is produced by sympa- 
thetic cooling with bosonic 23 Na atoms in a magnetic 
trap as described in [24]. The 6 Li atoms are then loaded 
into a single-beam optical dipole trap and transferred 
into the lowest hyperfme state \F = 1/2; tuf = 1/2) 
by a Landau-Zener radio-frequency (RF) sweep. Ad- 
ditional axial confinement is provided by magnetic 
field curvature. An equal mixture of |1) and |2) spin 
states (corresponding to the \F = 1/2; mp = 1/2) and 
\F = 1/2; = —1/2) states at low magnetic field) is 
prepared by a RF sweep at 300 G, followed by 500 ms 
wait time for decoherence and evaporative cooling in the 
optical trap. The Feshbach resonance at 834 G [1] is 
used to tune the repulsive interactions between |1) and 
|2). We increase the magnetic field in 200 ms to 528 
G, where the scattering length is zero and our Fermi 
gas is non-interacting. The final trap has a depth of 
4.1 /jK and frequencies of ui x = oj y = 390 Hz and uj z — 
34.7 Hz. The number of atoms per spin state is 8 x 10 5 , 
which corresponds to a Fermi temperature Tp of 1.4 
/xK. The temperature of the atoms is 0.3 Tf at this 
point. For loss rate measurements, the magnetic field is 
quickly ramped to the target field. For compressibility 
measurements, the field is ramped up over 50 ms and 
held for 30 ms to ensure thermal equilibrium before 
imaging. The molecular fraction in the density profile 
is determined by dissociating molecules with a magnetic 
field jump to 20 G above the Feshbach resonance, and 
comparing with the atom number after jumping to 528 
G where the cross section for imaging molecules vanishes 
(see Ref. gH). 



III. LOSS RATE 
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FIG. 1. Characterizing atomic loss for increasing repulsive 
interactions, (a) Remaining fraction of atom number imme- 
diately after the fastest possible ramp to the target field, (b) 
Atom loss rates at the target fields. Dashed line shows the 
estimated maximum tolerable loss rate of 0.02 /ms. (c) Molec- 
ular fraction after the 50 ms ramp and 30 ms wait time, cor- 
responding to when we measure equilibrium atomic density 
profiles. 



To obtain thermodynamic parameters from atomic 
density profiles requires equilibration. One time scale for 
equilibration is set by the longest trap period, which is 30 
ms for the axial direction. Ref. g| studied the distortions 
of profiles of bosonic 7 Li due to non-adiabatic ramping of 
the scattering length. The authors found that ramping 
times longer than rs 6tt/w (which is 90 ms for our pa- 
rameters) led to only small deviations from equilibrium 
profiles. Assuming that losses sharply increase during 
the last 5 ms of the ramp towards stronger interactions 
and limiting tolerable losses to 10 % leads to an estimate 
of a maximum tolerable loss rate of about 0.02/ms. The 
fastest relaxation time for excitations created during a 



non-adiabatic ramp is 1/ui (the damping time for a har- 
monic oscillator at critical damping). Allowing ~ 10 % 
loss during this time, leads to an identical estimate for 
the tolerable loss rate of ~ 0.02/ms. 

We measured loss rate as a function of interaction 
strength kp°a. Note that kp° is the Fermi wave vector 
of the zero temperature noninteracting gas calculated at 
the trap center using total atom number. The real kp 
is somewhat smaller because repulsive interactions and 
non-zero temperature lower the density. 

First, in Fig. QJa) we measured the number of atoms 
right before and after the fastest possible ramp (limited 
to 7 ms by eddy currents) to the target magnetic field. 
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During the ramp to the target magnetic field of kp°a 
0.8, ~ 35 % of the sample is lost. Measuring the loss rate 
at higher kp°a is difficult because most of the sample is 
lost even before reaching the target fields. The loss rate 
was determined by monitoring the drop in the number 
of atoms immediately after the field ramp. The results 
in Fig. [IJb) show that the measured loss rate reaches the 
maximum tolerable value of 0.02/ms at kp°a ss 0.35, lim- 
iting our measurements of equilibrium density profiles to 
smaller values of kp°a. Furthermore, at the same values 
of kF°a, the molecular fraction when we measure equilib- 
rium density profiles exceeded 10 % (see Fig.[IJc)). As we 
discuss in Sec.|lV]the presence of a molecular component 
affects the compressibility measurement. 



IV. COMPRESSIBILITY MEASUREMENT 

A system is characterized by its equation of state. The 
equation of state can be expressed in different forms in- 
volving various thermodynamic variables including den- 
sity, energy, pressure, temperature and entropy. For cold 
atom experiments, density, chemical potential (through 
the trapping potential) and temperature are directly 
accessible to measurement. In the weakly interacting 
regime, the interaction manifests itself as a perturbative 
term in the equation of state. Here we measure this per- 
turbative interaction effect by measuring the derivative 
of density with chemical potential, the isothermal com- 
pressibility of the gas. We prepare the system at the 
lowest temperature, but due to heating of the cloud by 
molecule formation, we have to apply a correction for the 
measured finite temperature. 

Using the experimental procedure discussed above we 
prepared equilibrated clouds at various magnetic fields 
and measured line-of-sight integrated profiles of column 
density using in situ phase-contrast imaging [26]. The 
signal-to-noise was improved by averaging the column 
density along equipotential lines (which are ellipses for 
the anisotropic harmonic oscillator potential). The av- 
eraging region was restricted to an axial sector of ±60° 
to avoid corrections due to transverse anharmonicities 
27J. Three-dimensional density profiles n(r) were recon- 
structed by applying the inverse Abel transformation to 
the column densities h(r) [26j . 

The isothermal compressibility is obtained from a spa- 
tial derivative of n(r), since in the local density approxi- 
mation, the local chemical potential is /i = \i Q — mui 2 x 2 /2, 
where /i Q is the global chemical potential and uj the trap 
frequency. The compressibility is defined as 
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We normalize the compressibility k by the compressibil- 
ity of an ideal gas at the given density and zero temper- 
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Here n is density, and m is the atomic mass. The normal- 
ized compressibility is obtained as the slope in a plot of 
n 2 / 3 vs. [i (Fig. 2). This plot is in essence the observed 
density profile plotted with the central region to the left 
and the spatial wings to the right. Experimentally, we 
find the slope to be constant over an extended range of 
the density profile. Compressibilities were determined 
from fits to the slope in the region of 90 to 50 % of the 
peak density. The region near the center of the cloud was 
excluded since the center is singular for the inverse Abel 
transformation leading to excessive noise. These com- 
pressibilities should be regarded as average values over 
the density range used in the fit. The uncertainty of fit- 
ting the slope to a single profile was 4.5 %. By averaging 
the profiles obtained from 20 images, the uncertainty was 
improved to 1.3 %. 

The normalized compressibility is a function of T/Tp 
and kpa. At finite temperature and scattering length 
a, T/Tp and kpa change across a single density profile 
because Tp and kp depend on the local density. However, 
this density dependence is small near the center of the 
cloud. Simulated density profiles showed that the average 
compressibility determined in the way described above 
agrees to within 0.6 % with the compressibility at T/Tp 
and kpa at the density in the center of the selected range. 
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FIG. 2. Determination of the compressibility of repulsively 
interacting Fermi gases. The compressibility is the slope of a 
graph showing the density to the power 2/3 versus the square 
of the distance from the center, (a) Single shot density profile, 
(b) Close-up of the region used for determining the compress- 
ibility for a plot showing the average of 20 density profiles. 



Compressibility decreases for stronger repulsion, but 
also for higher temperature. To identify the effect of 
repulsive interaction requires a careful consideration of 
finite temperature effects. First, the temperature of the 
cloud had to be accurately determined. This can be done 
without any special assumptions by fitting the wings of 
the cloud using a virial expansion |28l|. by thermome- 
try with another co-trapped atom [14j , or for population 
imbalanced clouds by fitting the wings of the majority 
component which is an ideal Fermi gas [29]. Here we 
chose to determine temperature using a virial expansion, 
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where X=JSK 



+ b 2 e 2 ^ + 0(e 3W ), 



(3) 



-, — ^ is the thermal de Broglie wavelength, 

b 2 is the virial coefficient, and is the fugacity. The 
virial coefficient for the Fermi gas with repulsive interac- 
tion is b 2 = — 2~ 5 / 2 — a/A [3(| ■ Pressure p was obtained 
from the doubly- integrated density profiles [f3[ . Temper- 
ature was determined in the wings of the profile where 
< —0.5. Here, the temperature measured with and 
without the interaction term —a/ A in b 2 differ by about 
3%. This suggests that higher-order corrections from in- 
teraction term will be negligible. Note that the virial 
expansion up to second order is valid to within 1% for 
the ideal gas at the density of the fitted wings. 

The low temperature normalized isothermal compress- 
ibility of a non-interacting Fermi gas is given by the Som- 
merfeld expansion [311 ] . 



«o,T =1 



T 



(4) 



where T is temperature and Tp is Fermi temperature. To 
add the effects of interactions, it is useful to work with 
the inverse normalized compressibility, 
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This is a derivative of the chemical potential, which has 
the following expansion in temperature and scattering 
length, 
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where C is constant, independent of density [321 ] . There- 
fore the inverse normalized compressibility has additive 
correction terms for temperature and interactions up to 
the second order of the interaction effect, 
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This equation defines Y(kpa), the interaction correc- 
tion to the inverse compressibility. This term is the 
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FIG. 3. Measured temperature, normalized compressibility, 
and interaction effect on compressibility at various interac- 
tion strengths, (a) The measured temperature as a function 
of interaction strength, (b) Solid circles show the measured 
normalized compressibility (k) and open squares show the 
calculated normalized compressibility at the measured tem- 
perature without interaction (k ,t). The difference between 
the two indicates the effect of interaction. Open triangles 
show the calculated normalized compressibility using the sec- 
ond order perturbation theory, which is consistent with our 
measured k. (c) The measured interaction correction to the 
inverse compressibility (solid circle) is compared to a linear fit 
(solid line), and the first (dashed line) and the second (dotted 
line) order perturbative result. 



derivative of the interaction term of the chemical poten- 
tial. In second order perturbation theory, one obtains 

vll \ 2; , 8(ll-2/n2)/ ; \ 2 

Y{k F a) = ±k F a + 1Sff 2 J {k F a,y. 

Fig. [3] shows the normalized compressibility, the tem- 
perature T/Tp, and the interaction correction to the in- 
verse compressibility. The temperature increases with 
k F a due to the increase in three-body recombination 
where the binding energy of the molecules (h 2 /ma 2 ) is 
transferred to the remaining atoms. The measured tem- 
perature is higher than that in previous experiments on 
the repulsive side [13, Hl|. This difference can be ex- 
plained by a smaller k F a since the increase in T/Tp is 
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approximately proportional to 1/kpa [21j |. 

We perform a linear fit of the interaction effect Y(kpd) 
versus k F a (constrained to pass through the origin) and 
obtain 0.680±0.147 for the slope, in agreement with the 
perturbative prediction of ^ = 0.637. This is the first 
observation of the mean-field term for repulsively inter- 
acting fermions in a thermodynamic quantity. The re- 
pulsive interaction has been seen as line shifts in RF 
spectroscopy experiments (which, in contrast to many 
thermodynamic quantities, are measured independently 
of the kinetic energy of the atoms) [33l. l34| . In principle, 
it is possible to obtain the mean-field term directly by 
fitting the density profiles with an extra mean field term. 
In such fits, we obtained clear evidence for such a term, 
but with low accuracy. It appears that the averaging over 
profiles for determining the compressibility(as in Fig. [2]) 
is superior. Figure [3] (c) shows the predicted effect of the 
second-order term on Y(kpa). With some improvements 
in signal-to-noise ratio, one should be able to observe this 
term which is analogous to the Lee-Huang- Yang correc- 
tion for bosons. 



V. DISPERSIVE EFFECT IN 
PHASE-CONTRAST IMAGING 

As mentioned in the introduction, phase-contrast 
imaging has several advantages over resonant absorp- 
tion imaging, and it has been app lied to many studies 
of cold Bose and Fermi gases [Is HE HH- Absorption 
imaging is usually done with absorptively dilute clouds, 
typically with 10 to 70 % absorption (or optical den- 
sities OD < 1). The standard assumption has been 
that dispersive imaging is quantitative when the phase 
shift 4> across the cloud is less than 7r/4. The normal- 
ized phase-contrast signal (for negligible absorption) is 
3 — 2v / 2cos(</) ± 7r/4), which is equal to 1 ± 2(f) for small 
phase shifts [23[ . The sign depends on laser detuning and 
the sign of the phase shift imparted by the phase plate. 

Here, we apply phase-contrast imaging for rather pre- 
cise quantitative studies of ultracold Fermi gases and 
found that even for small phase shifts systematic dis- 
persive distortions of the image cannot be neglected. 
Phase-contrast imaging relies on column density depen- 
dent phase shifts. However, if the object is not thin, but 
extended along the line of sight, some lensing will affect 
the images. These distortions should vary inversely pro- 
portional to the probe light detuning and become negli- 
gible for far detuning. 

We investigated positive and negative detuning (r5) 
of 40, 80, and 120 MHz (see Fig. IJ). The normalized 
phase contrast image had a maximum signal of 0.35/1.85, 
0.55/1.6, and 0.7/1.4 for the three positive/negative de- 
tuning. The lensing effect is opposite (focusing vs. defo- 
cusing) for positive and negative detuning, and can there- 
fore be identified by comparing profiles obtained with 
positive and negative detuning. Fig. [5] shows that at 
40 MHz, the two profiles show a visible difference, but 




FIG. 4. Phase-contrast imaging of a balanced spin mixture 
in states |1) and |2). White images (phase shift <j> >0) were 
obtained for a probe beam red-detuned from the |2) — s- |e) 
transition, corresponding to 5 < 0. Black images (<j) <0) were 
obtained for a probe beam blue-detuned from the |1) — > |e) 
transition, corresponding to 8 < 0. 
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FIG. 5. Comparison of column density profiles obtained from 
positive and negative detuning of (a) 40 and (b) 80 MHz, 
respectively, (c) The measured compressibility (at fc^a = 
0) from positive (open triangle) and negative (open square) 
detuned profiles and their averages (solid circle) are shown. 
The average value stays constant above 80 MHz detuning. 



for profiles at 80 MHz and 120 MHz, the differences are 
smaller than the noise level. 

However, the compressibility is determined by the 
slope of the profiles and very sensitive to small distortions 
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even if they are not perceptible in the profiles. Fig. [5jc) 
shows that even at 120 MHz detuning, the compress- 
ibilities obtained from profiles with the two signs of the 
detuning differ by about 10 %. Since further detuning 
would have resulted in a smaller signal we evaluated the 
average value of the compressibility for positive and neg- 
ative. When the dispersive distortions are small, the ef- 
fect on the compressibility should be a first-order effect 
in the phase shift and cancel for the average. Indeed, the 
average value stays constant above 80MHz detuning. Our 
conclusion is that for reasonable signal levels (i.e. 50 % 
of the baseline set by the probe light) dispersive effects 
are relevant for quantitative studies, but can be elimi- 
nated by performing averages over positive and negative 
detunings. 



VI. DISCUSSION 

We address now to what extent a small molecular frac- 
tion contributes to the observed density profiles. The 
presence of molecules is unavoidable since they form dur- 
ing the ramping and equilibration time. At the highest 
magnetic field used in the experiment, 679 G, the molec- 
ular fraction was determined to be ~ 10 % (Fig. 1). Ref. 
[37| reported that molecules at 650G showed an absorp- 
tion cross section of about half the value of the atoms for 
probe light at the atomic resonance. For phase-contrast 
imaging with large detuning molecules and atoms should 
contribute equally to the signal. 

We performed simulations to address how the presence 
of molecules would affect the compressibility measure- 
ments. We considered as possible scenarios (i) that the 
molecular fraction is constant throughout the cloud, (ii) 
that the molecular fraction is proportional to the loss rate 
(n 8 / 3 ), and (iii) that the molecular fraction is well equi- 
librated at the same temperature as atoms. The atomic 
profile is then the difference of the measured density pro- 
file minus the simulated molecular density distribution. 
Scenario (iii) is ruled out since it would result in a rather 
sharp peak in the density profile which was not observed. 
The first two scenarios with a 10 % molecular fraction re- 
sulted in a value for the normalized compressibility which 
was increased by 3.3 % and 4.4 % respectively. This 
shows that for our largest value of kpa the presence of 
molecules starts to become a systematic effect. In addi- 
tion to the contribution to the density profiles, molecules 
can affect the atomic density distribution through atom- 
molecule interactions. 

Our work shows that the interaction effect on the com- 
pressibility at the maximum possible values of kpa is 
about 15 %. We could identify this effect only by careful 
thermometry (to distinguish it from thermal effects) and 
by correcting small dispersive distortions of the cloud. 
It is desirable to study fermions for stronger repulsive 
interactions where stronger and non-perturbative effects 
are predicted. The maximum possible kpa value for 
obtaining equilibrium density profiles is determined by 



the loss rate which is proportional to n 2 a 6 max(T, Tp) = 
(k F a) 6 n 2 / 3 0, HI. Therefore, the maximum possible 
kpa for a given loss rate is proportional to n~ x / 9 and 
stronger interaction effects can be seen at lower density. 
This should be accomplished by reducing the radial con- 
finement and not the axial confinement which determines 
the equilibration time. However, the weak density depen- 
dence will allow only modest increases in kpa. A recent 
experiment used density ten times smaller than ours [2l| 
and reported ramping from kp°a = to kp°a = 0.35 in 
500 ms losing only 5 % atoms. Assuming the loss hap- 
pened during the last 50 ms, we can roughly estimate a 
loss rate of ~ 0.001 ms -1 at kp°a = 0.3 which is lower 
than our measurement, consistent with the lower density. 

Longer lifetimes for a given kpa should be realized 
using narrow instead of broad Feshbach resonances. 
For narrow Feshbach resonances the low-lying molecular 
state has a dominant closed-channel character. There- 
fore, three body recombination of atoms (which are in 
the open channel) has a smaller overlap to the molecular 
state and therefore a reduced loss rate. Recent experi- 
ments using RF spectroscopy [Z(| El[ confirm this. How- 
ever, for such narrow resonances the zero-range approx- 
imation is no longer valid, the interaction is no longer 
described by the scattering length alone and becomes 
(through an effective range parameter) momentum de- 
pendent. As a result, the narrow Feshbach resonances 
realize a different Hamiltonian. 

In conclusion, in this paper we have addressed the 
question to what extent Fermi gases with strong inter- 
actions can be studied by observing equilibrium density 
profiles. The range of sufficiently long metastability to 
reach equilibrium is limited to values of/ci?a<0.25. In 
this range, interaction effects are comparable to thermal 
effects, but we were able to observe how interactions re- 
duce the compressibility and obtained quantitative agree- 
ment with the first-order mean field term. An observa- 
tion of the second order Lee-Huang- Yang correction is in 
experimental reach. 

If experiments can be performed at stronger interac- 
tions, a natural extension of our work would be a mea- 
surement of the spin susceptibility using population im- 
balanced Fermi systems. This was performed recently for 
fermions with attractive interactions [15| . Such measure- 
ments could address the possible existence of a ferromag- 
netic transition in a repulsive Fermi gas [20| for which the 
spin susceptibility would diverge at the phase transition 
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